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FOREWORD 


The  1995  version  of  the  Aeroprediction  Code  (AP95)  extended  the  code’s  angle-of-attack 
(AOA)  capability  to  90  deg  through  empirical  additions  to  normal  force  of  various  missile 
components  based  on  several  large  wind  tunnel  data  bases.  However,  very  little  emphasis  was  placed 
on  improving  the  axial  force  coefficient  at  high  AOA.  The  justification  for  this  was  that  normal  force 
was  by  far  the  dominant  term  and  axial  force  only  a  minor  contributor  to  lift  and  drag.  While  this 
statement  is  still  true,  it  is  less  true  at  subsonic  speeds  and  therefore  can  result  in  appreciable  errors 
in  axial  force  and  drag  prediction  at  high  AOA.  As  a  result,  new  technology  has  been  developed  to 
remedy  this  problem.  This  new  technology  is  discussed  in  this  report. 

The  work  described  in  this  report  was  supported  through  the  Office  of  Naval  Research 
(Mr.  Dave  Siegel)  by  the  following  programs:  the  Air  Launched  Weapons  Program  managed  at  the 
Naval  Air  Warfare  Center,  China  Lake,  CA,  by  Mr.  Tom  Loftus  and  Dr.  Craig  Porter,  and  the 
Surface  Weapons  Systems  Technology  Program  managed  at  the  Naval  Surface  Warfare  Center, 
Dahlgren  Division  (NSWCDD)  by  Mr.  Robin  Staton  and  Mr.  Gil  Graff  Also,  some  support  was 
provided  by  the  Marine  Corps  Weaponry  Technology  Program  managed  at  NSWCDD  by  Mr.  Bob 
Stiegler.  The  authors  express  appreciation  for  support  received  in  this  work. 


Approved  by: 


DAVID  S.  MALYEVAC,  Deputy  Head 
Weapons  Systems  Department 


iii/iv 


NSWCDD/TR-96/240 


CONTENTS 

Section  Page 


1.0  INTRODUCTION .  1 

2.0  ANALYSIS .  2 

2.1  REVIEW  OF  C^  METHODOLOGY  IN  THE  AP95  .  2 

2.2  Ca  PHYSICAL  PHENOMENA  AT  AOA .  3 

2.3  IMPROVED  SEMIEMPIRICAL  AOA  THEORETICAL  MODEL .  4 

3.0  RESULTS  AND  DISCUSSION .  15 

4.0  SUMMARY  AND  RECOMMENDATIONS .  39 

5.0  REFERENCES  .  42 

6.0  SYMBOLS  AND  DEFINITIONS .  44 

DISTRIBUTION .  (1) 

ILLUSTRATIONS 

Figures  Page 

1  AXIAL  FORCE  AOA  VARIATION  PARAMETERS  FOR  BODY  ALONE  . .  6 

2  AXIAL  FORCE  AOA  VARIATION  PARAMETERS  FOR  A  BODY-TAEL 

CONHGURATION  .  8 

3  AXIAL  FORCE  AOA  VARIATION  PARAMETERS  FOR  A  WING-BODY- 

TAIL  CONHGURATION  .  10 

4  f(M,aw)  VALUES  FOR  a  AND  6  OF  OPPOSITE  SIGNS  .  12 

5A  BODY  ALONE  AND  ONE  BODY-TAIL  CONHGURATION  OF  NASA 

TRI-SERVICE  DATA  BASE  .  16 

5B  COMPARISON  OF  AXIAL  FORCE  COEFFICIENTS  OF  THEORY  AND 

EXPERIMENT  FOR  NASA  BODY  ALONE  CONHGURATION 
OF  FIGURE  5A .  17 


V 


r 


NSWCDD/rR-96/240 
ILLUSTRATIONS  (Continued) 

Figures  Page 

5C  COMPARISON  OF  AXIAL  FORCE  COEFFICIENTS  OF  THEORY  AND 

EXPERIMENT  FOR  NASA  BODY-TAIL  CONFIGURATION 

OFHGURESA . .  20 

6A  BODY  ALONE  AND  ONE  BODY-TAIL  CONHGURATION  OF 

BAKER  DATABASE'^  .  22 

6B  COMPARISON  OF  AXIAL  FORCE  COEFFICIENTS  OF  THEORY  AND 

EXPERIMENT'^  FOR  BODY  ALONE  OF  FIGURE  6A .  23 

6C  COMPARISON  OF  AXIAL  FORCE  COEFHCIENTS  OF  THEORY  AND 

EXPERIMENT'^  FOR  BODY-TAIL  CONFIGURATION  OF 

FIGURE  6A .  25 

7A  CANARD-BODY-TAIL  CONnGURATION  WITH  HEMISPHERICAL 

NOSE'’  .  27 

7B  AXIAL  FORCE  VERSUS  AOA  FOR  THE  CONFIGURATION  OF 

HGURE  7A .  28 

8A  CANARD-CONTROLLED  MISSILE  CONHGURATION'*  (WIND 

TUNNEL  MODEL  1/3  SCALE) .  30 

8B  COMPARISON  OF  THEORETICAL  METHODS  FOR  AXIAL  FORCE 

PREDICTION  TO  EXPERIMENT  ON  HGURE  8A  CONFIGURATION 

($  =  0  deg,  M  =  0.2,  Rn=  1.42  X  10-Vft)  .  31 

8C  COMPARISON  OF  THEORETICAL  METHODS  FOR  AXIAL  FORCE 

PREDICTION  TO  EXPERIMENT  ON  FIGURE  8A  CONFIGURATION 

(O  =  45deg,M  =  0.2,RN=  1.42X  lO'Vft)  .  32 

9A  AIR-TO-AIR  MISSILE  CONHGURATION  USED  IN  VALIDATION’® .  34 

9B  COMPARISON  OF  THEORY  AND  EXPERIMENT  FOR  CONFIGURATION 

OFHGURE9A(<I>  =  Odeg,6  =  Odeg)  .  35 

9C  COMPARISON  OF  THEORY  AND  EXPERIMENT  FOR  CONFIGURATION 

OFFIGURE9A(€>  =  45deg,6  =  Odeg)  .  36 

9D  COMPARISON  OF  THEORY  AND  EXPERIMENT  FOR  CONHGURATION 

OF  HGURE  9A  (<I>  =  0  deg) .  37 

9E  COMPARISON  OF  THEORY  AND  EXPERIMENT  FOR  CONHGURATION 

OF  FIGURE  9A  (O  =  45  deg) . . .  38 

lOA  COMPARISON  OF  AXIAL  FORCE  COEFHCIENTS  OF  A  TAIL 

CONTROLLED  WING-BODY-TAIL  CONFIGURATION  (0  =  0  deg) ...  40 

lOB  COMPARISON  OF  AXIAL  FORCE  COEFHCIENTS  OF  A  TAIL 

CONTROLLED  WING-BODY-TAIL  CONFIGURATION  (O  =  45  deg) . .  41 


VI 


NSWCDD/TR-96/240 


TABLES 

Tables  Page 


1  VALUES  OF  THE  AXIAL  FORCE  AOA  PARAMETERS  FOR  A 

BODY  ALONE  CONHGURATION .  6 

2  VALUES  OF  THE  AXIAL  FORCE  AOA  PARAMETERS  FOR  A 

BODY-TAIL  CONFIGURATION  .  9 

3  VALUES  OF  THE  AXIAL  FORCE  AOA  PARAMETERS  FOR  A 

WING-BODY-TAIL  CONFIGURATION .  11 

4A  f(M,aJAT$  =  Odeg .  13 

4B  f(M,a  J  AT  O  =  45  deg .  14 


vii 


NSWCDD/TR-96/240 


1.0  INTRODUCTION 


The  latest  version  of  the  Aeroprediction  Code  (AP95)*  available  to  the  public  places  a  lot  of 
emphasis  on  zero  lift  drag  or  axial  force  coefficient  at  zero  angle  of  attack  (AO A)  with  only  minor 
attention  devoted  to  the  change  in  C^  with  AOA.  One  component  of  C^,  the  base  axial  force 
coefficient  change  with  AOA,  is  treated  in  a  fairly  comprehensive  manner.  The  primary  rationale  for 
not  devoting  more  effort  to  calculating  the  change  in  total  axial  force  with  AOA  is  the  fact  that  for 
most  conditions,  normal  force  is  the  dominant  term  in  lift,  drag,  and  lift/drag  ratio  at  high  AOA,  and 
Ca  plays  a  secondary  role.  However  at  subsonic  Mach  numbers,  C^  can  actually  go  negative  at  high 
AOA,  causing  more  error  in  Cl,  Cp,  and  I7D  than  desired.  Also,  at  very  high  Mach  numbers 
(M  >  3.0),  Ca  tends  to  increase  with  AOA  much  more  than  the  current  AP95  predicts.  Finally,  the 
change  in  Ca  with  control  deflection,  while  predicted  reasonably  well  by  the  AP95,  could  possibly  be 
improved  somewhat,  also  at  subsonic  and  high  supersonic  conditions.  It  is  therefore  the  intent  of  this 
report  to  discuss  the  physics  of  the  changes  in  axial  force  with  AOA,  develop  a  semiempirical 
mathematical  model  to  account  for  the  physics,  and  compare  the  improved  axial  force  model  at  AOA 
with  the  AP95  and  experimental  data. 

In  researching  the  literature  for  methods  to  account  for  axial  force  changes  with  AOA,  three 
were  found.  The  first  of  these  was  by  Jorgensen^  where  he  approximated  Ca  by: 

Ca  “  Ca^  cos^a  ;  0  <  a  <  90°  (1) 


This  method  worked  reasonably  well  for  subsonic  Mach  numbers  and  AO  As  less  than  about  30  deg. 
References  3  and  4  improved  upon  the  Jorgensen  method  by  assuming 

Ca  =  C^^  .  f(M,a)  (2) 


f(M,a)^  was  approximated  by  a  fourth  order  polynomial  in  AOA  for  subsonic  and  transonic  Mach 
numbers  and  was  assumed  to  be  zero"  for  supersonic  Mach  numbers.  This  was  an  improvement  on 
the  Jorgensen  methodology  at  all  Mach  numbers  and  particularly  for  AOAs  greater  than  about 
30  deg.  However,  since  f(M,a)  =  0  at  M>  1.5,  the  method  did  not  pick  up  the  increase  in  Ca  due  to 
compressibility  effects  at  higher  M.  Reference  5  then  improved  upon  Reference  4  by  assuming  f(M,a) 
was  second  order  in  AOA  at  supersonic  Mach  numbers  for  bodies  alone  and  third  order  for  wing- 
bodies.  Unlike  Reference  1,  which  gives  consistently  accurate  results  for  Ca  ,  using  second  order 
methods,  References  2  through  5  all  use  basically  empirical  methods  for  computing  .  This  means 

the  empirical  methods  are  limited  to  the  data  bases  upon  which  they  are  based,  not  only  at  high  AOA, 
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but  low  AOA  as  well.  Furthermore,  neither  References  4  or  5  treated  wing-body-tail  configurations 
separately  from  wing-body  cases. 

While  References  4  and  5  give  acceptable  accuracy  for  f(M,a)  at  AOAs  to  90  deg  for  body  alone 
and  body-tail  cases,  they  use  different  prediction  methods  in  various  Mach  number  regions.  It  is  the 
goal  of  the  present  work  to  derive  a  single  new  function,  f(M,a),  that  is  as  accurate  or  more  accurate 
than  either  the  Reference  4  or  5  methods,  yet  is  applicable  over  the  entire  Mach  number  and  AOA 
range  of  interest  and  is  applicable  to  body  alone,  wing-body,  or  wing-body-tail  cases.  In  addition, 
neither  Reference  2, 4,  or  5  treated  empirical  corrections  to  axial  force  on  missile  configurations  with 
combined  AOA  and  control  deflection.  This  area  will  also  be  investigated  for  possible  improvement. 


2.0  ANALYSIS 


2.1  REVIEW  OF  C,  METHODOLOGY  IN  THE  AP95 

Equation  (2)  will  be  assumed  to  apply  for  both  body  alone,  and  configurations  with  lifting 
surfaces  as  well.  The  AP95  methodology  for  of  Equation  (2)  is  based  on  computing 
independently  the  wave  or  pressure  axial  force,  skin-friction  and  base  term.^  The  wave  drag  term  is 
computed  at  low  supersonic  Mach  numbers  by  combining  modified  Newtonian  theory  (MNT)  with 
Second-Order  Van  Dykes  hybrid  theory,^  at  moderate  supersonic  Mach  numbers  by  improvements 
to  second-order,  shock  expansion  theory  combined  with  MNT,’  and  at  Hypersonic  Mach  numbers 
by  further  refinements  in  shock  expansion  theoiy  and  MNT  extended  to  real  gases.*  The  combination 
of  the  above  new  methods  developed  during  the  course  of  the  AP95  development  has  allowed  zero 
lift  drag  to  be  computed  with  average  accuracies  generally  of  ±10  percent  for  a  broad  range  of 
configurations  and  for  Mach  numbers  ^1.2.  Configurations  can  have  truncated,  spherically  blunt  or 
sharp  noses  with  boattails  or  flares.  Additionally,  up  to  two  discontinuities  in  surface  slope  are 
allowed  along  the  nose  of  the  body. 

At  subsonic  Mach  numbers,  the  pressure  component  of  axial  force  is  assumed  to  be  zero  except 
in  cases  where  a  large  cone  half  angle  exists  on  the  nose,  creating  a  separation  bubble  in  the  vicinity 
of  the  nose  or  on  truncated  noses.  The  viscous  separation  term  is  estimated  empirically.**  At 
transonic  Mach  numbers,  the  wave  drag  is  estimated  analytically  by  either  the  method  of  Wu  and 
Aoyoma’  or  using  an  unsteady  solution  of  the  Euler  equations.*®  These  estimates  of  wave  drag  were 
then  included  in  table  lookup  form  into  the  AP95®’**  as  functions  of  nose  length,  bluntness  and  Mach 
number.  Boattail  effects  were  computed  by  the  Reference  9  methodology. 

The  skin-friction  drag  was  computed  by  the  Van  Driest  n  method”  for  the  turbulent  portion  of 
the  boundary  layer.  The  laminar  portion  is  based  on  the  incompressible  flow  over  a  flat  plate  and 
modified  for  compressibility  in  Reference  1 .  The  AP95  allows  various  options  for  surface  conditions 
including  typical  flight,  wind  tunnel  model  with  and  without  a  boundary  layer  trip,  and  all  laminar 
flow.  These  options  affect  the  transition  Reynolds  number  where  the  flow  becomes  turbulent. 
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The  base  drag  is  computed  empirically  as  a  function  of  AO  A,  Mach  number,  fin  location,  fin 
thickness  effects,  fin  control  deflection,  and  boattail  and  for  both  power  on  or  off  conditions.  The 
methodology  is  given  in  Reference  1 3  and  was  upgraded  from  the  original  AP72®  by  new  wind  tunnel 
data  described  in  Reference  13. 

Fin  effects  on  are  computed  in  a  similar  way  to  the  body  alone.  There  are  a  couple  of 
differences  however,  for  blunt  leading  edge  wings,  the  MNT  is  combined  with  first  order  versus 
second  order  perturbation  theory.  This  can  be  done  for  two  reasons.  First,  the  wing  drag  is  generally 
much  smaller  than  the  body  drag.  Secondly,  the  wings  generally  are  thin,  thus  having  frontal  areas 
with  smaller  slopes  than  the  body.  This  allows  first  order  perturbation  theory  to  yield  acceptable 
wave  drag  on  wings.  The  second  difference  between  the  of  wings  and  bodies  is  the  wing  trailing 

edge  drag  versus  body  base  drag.  The  body  base  drag  uses  a  three-dimensional  base  pressure 
coefficient  whereas  a  two-dimensional  base  pressure  coefficient  is  used  for  a  blunt  wing  trailing  edge. 

As  is  apparent  from  the  discussion  on  ,  the  AP95  has  placed  a  major  emphasis  on  accurate 
values  of  from  its  inception  in  1972.  The  reason  for  this  goes  back  to  the  fact  that  in  1972,  there 
was  no  accurate  way  to  calculate  drag  of  spin-stabilized  projectiles  for  range  estimation.  To  put  the 
level  of  emphasis  on  in  context  of  the  overall  AP95,  consider  the  following.  The  AP95  contains 
135  primary  subroutines  consisting  of  over  16000  lines  of  code.  It  is  estimated  that  about  half  the 
code  is  for  computing  accurate  values  of  .  One  simply  cannot  get  accurate  values  of  for  a 
broad  class  of  body  shapes  by  empirical*^ or  first  order  perturbation  methods.  Seconcf-order 
perturbation  methods,  in  combination  with  modified  Newtonian  theory,  and  judicious  selection  of  a 
match  point  between  Newtonian  and  the  second-order  method  is  the  lowest  order  analytical  method 
acceptable  for  projectile  and  missile  range  estimation  in  general. 


2.2  Ca  physical  PHENOMENA  AT  AOA 

Before  developing  a  mathematical  model  to  address  changes  in  axial  force  coefficient  as  a 
function  of  AOA,  it  is  appropriate  to  discuss  the  physics  of  the  flow  that  causes  these  changes.  To 
visualize  the  changes,  the  axial  force  will  be  once  again  broken  down  into  its  components  due  to 
pressure,  skin-friction  and  base  drag. 

For  subsonic  Mach  numbers  and  spherically  blunt  or  ogive  shaped  bodies,  the  pressure  drag  at 
low  to  moderate  AOA  is  zero.  However,  as  AOA  is  increased  or  for  bodies  with  truncated  noses  or 
large  cone  half  angles,  the  flow  forms  a  separation  bubble  in  the  vicinity  of  the  nose  region.  This 
separation  bubble  has  a  negative  pressure  coefficient,  which  means  the  axial  force  decreases  with 
AOA.  On  the  other  hand,  at  high  Mach  number,  the  pressure  coefficient  on  the  windward  side  of  the 
body  is  a  function  of  the  sine  squared  of  the  angle  between  a  tangent  to  the  body  surface  and  the 
velocity  vector.  This  means  the  increase  in  axial  force  at  high  Mach  number  with  AOA  is  positive 
due  to  compressibility  effects  of  the  air.  On  the  leeward  surface  of  the  body,  the  pressure  coefficient 
approaches  zero  at  high  Mach  number  so  it  has  little  effect  on  the  axial  force  coefficient.  None  of 
the  theoretical  methods  being  used  for  computing  the  pressure  component  of  the  axial  force 
coefficient  accurately  account  for  these  changes  above  AOA  of  about  10  deg. 
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The  base  axial  force  changes  with  AOA  are  positive  at  all  Mach  numbers’^  but  are  the  largest  at 
transonic  Mach  numbers.  At  very  high  Mach  numbers  where  the  base  term  approaches  zero,  the 
change  with  AOA  is  also  zero.  At  moderate  AOA,  the  increase  in  base  axial  force  reaches  a 
maximum  and  then  begins  decreasing  with  additional  increases  in  AOA.  When  fins  are  placed  on  the 
body  in  the  vicinity  of  the  base,  the  increase  in  base  axial  force  with  AOA  is  mitigated  substantially. 
It  is  not  clear  the  phenomena  that  cause  these  effects  with  AOA.  It  is  suspected  that  for  the  body 
alone,  the  initial  AOA  increases  cause  a  stronger  separation  at  the  shoulder  due  to  a  larger  turn  angle 
at  the  base.  However,  when  fins  are  placed  in  the  vicinity  of  the  base,  it  is  suspected  that  some  of  the 
high  dynamic  pressure  fluid  is  channeled  into  the  base  region,  thus  mitigating  the  negative  base 
pressure  coefficient  somewhat. 

The  skin-friction  coefficient  axial  force  change  with  AOA  is  even  harder  to  estimate  or  to 
explain,  due  to  the  lack  of  CFD  or  wind  tunnel  data  which  attempt  to  measure  this  quantity  at 
moderate  to  high  AOA.  At  lower  Mach  number  where  strong  separation  exists  in  the  leeward  plane 
area  along  with  a  fairly  large  reverse  flow  region,  it  is  suspected  that  the  skin  friction  decreases  with 
AOA.  As  Mach  number  increases,  the  reverse  flow  region  decreases  and  the  change  in  skin  friction 
axial  force  with  AOA  is  probably  nearly  constant  or  maybe  slightly  positive  for  low  to  moderate 
AO  As. 

Another  physical  phenomena  which  causes  changes  in  (at  moderate  and  higher  supersonic 
Mach  numbers)  with  AOA  is  internal  shock  interactions.  These  phenomena  are  primarily  associated 
with  bodies  with  lifting  surfaces  and  particularly  those  with  more  than  one  set.  The  body  alone 
generally  will  have  only  a  bow  shock  and  internal  expansion  waves  unless  there  is  a  surface 
discontinuity  present.  A  body-tail  configuration  will  have  shocks  coming  off  the  tail  surfaces  which 
interact  with  the  bow  shock.  At  Mach  numbers  and  AOAs  high  enough,  these  interactions  can  have 
an  impact  on  the  axial  force.  Finally,  if  two  or  more  sets  of  lifting  surfaces  are  present,  even  stronger 
shock  interactions  occur  between  not  only  the  bow  and  wing  shock  but  wing  to  tail  shocks  as  well. 
The  situation  is  further  complicated  if  a  control  surface  is  deflected.  The  internal  shock  interaction 
effects  can  affect  the  other  static  aerodynamics  as  well  as  axial  force  when  this  is  the  case. 


2.3  IMPROVED  SEMIEMPIRICAL  AOA  THEORETICAL  MODEL 

As  stated  in  the  Introduction  section  of  the  report,  the  goal  is  to  develop  a  single  method  to 
predict  f(M,a)  of  Equation  (2)  as  accurate  as  either  of  the  methods  in  References  4  and  5,  but  which 
applies  for  all  Mach  numbers  for  bodies  with  and  without  lifting  surfaces,  and  for  AOA  to  90  deg. 
To  this  end,  we  assume 


f(M,a)  =  Aa  +  Ba^  +  Ca^  +  Da"*  (3) 

To  evaluate  the  coefficients  A,  B,  C,  and  D  require  four  independent  conditions.  These  conditions 
are: 


_af 

da 


=  f  (M,0) 


la=0 
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ii)-iv)  The  value  of  f(M,a)  at  a  =  30  deg,  60  deg  and  90  deg  respectively. 

Using  these  four  conditions,  and  putting  a  in  radians  versus  degrees,  four  equations  are  arrived  at  to 
solve  simultaneously.  These  equations  are: 


A 

.524A 

=  f(M,0) 

+  .274B 

+  .144C 

+  .076D 

=  f(M,30) 

1.047A 

+  1.096B 

+  1.148C 

+  1.202D 

=  f(M,60) 

1.571  A 

+  2.467B 

+  3.875C 

+  6.059D 

=  f(M,90) 

Simultaneous  solution  of  the  variables  A,  B,  C,  and  D  from  Equation  (4)  gives 
A=  f(M,0) 

B=  -3.509  f(M,0)  +11.005  f(M,30)  -2.757  f(M,60)  +0.41  f(M,90) 

C=  3.675  f(M,0)  -17.591  f(M,30)  +7.041  f(M,60)  -1.179  f(M,90) 

D-  -1.181  f(M,0)  +6.771  f(M,30)  -3.381  f(M,60)  +.752  f(M,90) 

Knowing  f'(M,0),  f(M,30),  f(M,60)  and  f(M,90)  in  conjunction  with  Equations  (5),  (3)  and  (2), 
improved  estimates  of  at  AOAs  up  to  90  deg  should  be  obtained.  In  reality,  if  the  body  is 
symmetric,  then  this  AOA  should  be  ±90  so  long  as  Equation  (3)  is  viewed  in  terms  of  absolute  values 
forAOA. 

The  question  that  must  be  addressed  is  how  to  determine  values  of  the  parameters  f'(M,0), 
f(M,30),  f(M,60)  and  f(M,90).  To  do  this,  two  large  wind  tunnel  data  bases  were  utilized.  These 
were  the  NASA  Tri-Service  Data  Base,*'*  which  contains  detailed  missile  component  data  to  AOA 
45  deg,  and  the  Baker^^  high  AOA  data  base,  which  includes  missile  component  data  to  180  deg 
AOA.  Reference  14  was  used  primarily  for  the  lower  AOA  information,  f'(M,0)  and  f(M,30),  and 
Reference  15  was  used  for  the  higher  AOA  variables,  f(M,60)  and  f(M,90).  The  reason  for  this  was 
that  Reference  14  data  were  available  in  tabular  form  making  the  estimates  of  the  parameters  more 
accurate  than  from  the  Reference  15  information,  which  was  available  only  in  graphs.  For  some  data 
points  in  the  transonic  flow  region,  the  parameter  f(M,30)  was  obtained  by  averaging  the  data  from 
both  References  14  and  15. 

Figure  1  gives  the  values  of  the  parameters  f'(M,0),  f(M,30),  f(M,60)  and  f(M,90)  for  the 
body  alone.  Reference  5  indicated  analysis  showed  the  first  parameter  (f'(M,0))  was  independent  of 
nose  length,  nose  shape,  and  total  body  length  and  dependent  only  on  Mach  number.  This  assumption 
will  also  be  made  for  the  other  three  parameters  as  well.  Comparison  of  the  new  method  on 
configurations  different  from  those  within  the  data  base  will  determine  the  validity  of  this  assumption. 
Table  1  gives  the  values  of  the  body  alone  parameters  of  Figure  1. 
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FIGURE  1 .  AXIAL  FORCE  AOA  VARIATION  PARAMETERS 
FOR  BODY  ALONE 


TABLE  1.  VALUES  OF  THE  AXIAL  FORCE  AOA  PARAMETERS  FOR 
A  BODY  ALONE  CONHGURATION 


M 

f(M,0)/rad 

f(M,30) 

f(M,60) 

f(M,90) 

0 

0.31 

-0.030 

-1.00 

0 

0.6 

0.31 

-0.030 

-1.00 

0 

0.8 

0.31 

-0.083 

-1.00 

0 

0.9 

0.31 

-0.095 

-0.81 

0 

1.0 

0.40 

-0.030 

-0.54 

0 

1.15 

0.67 

0.060 

-0.36 

0 

1.3 

0.65 

0.090 

-0.27 

0 

1.5 

0.63 

0.120 

-0.04 

0 

2.0 

0.49 

0.142 

0.19 

0.045 

2.5 

0.29 

0.147 

0.24 

0.090 

3.0 

0.20 

0.150 

0.30 

0.090 

4.5 

0.18 

0.173 

0.38 

0.090 

1 

2:6.0 

0.16 

0.210 

0.40 

0.090 
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Figure  2  and  Table  2  give  the  values  of  the  same  parameters  for  a  configuration  with  one  set 
of  lifting  surfaces.  These  parameters  were  derived  for  mostly  body-tail  configurations.  In  comparing 
values  of  the  parameters  in  Figure  2  with  those  of  Figure  1 ,  it  is  seen  that  similar  values  exist  for  each 
parameter  but  they  are  slightly  different.  The  Reference  4  method  shows  no  difference  between 
configurations  that  are  bodies  alone  versus  those  with  lifting  surfaces,  and  Reference  5  shows  a 
difference  only  as  AOA  approaches  90  deg.  Thus  the  present  method  will  also  give  slightly  different 
results  for  body  alone  and  wing-body  or  body-tail  missile  configurations  than  those  of  References  4 
and  5  due  to  differences  in  the  values  of  the  parameters  used  in  defining  the  nonlinearity  with  AOA. 
This  is  in  addition  to  differences  arising  from  the  use  of  a  single  equation  for  all  cases  in  the  present 
approach  versus  three  different  equations  used  in  References  4  and  5. 

Unfortunately,  neither  of  the  data  bases  given  in  References  14  or  15  tested  configurations 
with  two  sets  of  lifting  surfaces.  While  the  coefficients  defined  by  Figure  2  and  T able  2  could  be  used 
as  an  approximation  for  configurations  with  two  sets  of  lifting  surfaces,  it  will  be  attempted  to 
improve  upon  these,  at  least  for  the  lower  AOAs  where  other  data  are  available.  To  this  end. 
References  17  through  20  were  utilized.  References  17  and  18  were  veiy  helpful  in  defining  f'(M,0) 
and  f(M,30)  for  Mach  numbers  0.8  to  4.6.  Above  4.6,  values  of  these  parameters  were  extrapolated 
as  was  done  in  Tables  1  and  2  as  well.  References  19  and  20  were  utilized  for  low  Mach  number 
values  of  f'(M,0),  f(M,30)  and  f(M,60).  f(M,60)  for  high  Mach  numbers  could  be  extrapolated 
reasonably  well  based  on  the  Reference  18  data  in  conjunction  with  Figure  2.  The  f(M,90)  data  are 
simply  a  best  guess  based  on  f(M,60)  values  for  wing-body-tail  cases  and  f(M,60)  and  f(M,90)  values 
of  Table  2. 

Values  of  the  four  parameters  of  Equation  (5)  are  given  in  Figure  3  and  Table  3  for 
configurations  having  two  sets  of  lifting  surfaces.  In  comparing  Figures  2  and  3,  it  is  seen  there  are 
some  similarities  but  also  some  differences.  At  low  Mach  numbers,  trends  of  all  four  parameters  are 
similar,  which  means  that  Table  2  could  be  used  successfully  for  configurations  with  more  than  one 
set  of  lifting  surfaces  and  still  get  reasonably  accurate  estimates  of  axial  force  change  with  AOA.  This 
is  in  fact  what  References  4  and  5  both  do.  However,  as  Mach  number  increases,  values  of  the 
parameters  f'(M,0)  and  f(M,30)  tend  to  be  higher  for  two  sets  of  lifting  surfaces  compared  to  one. 
That  is,  Ca  increases  faster  and  reaches  a  higher  peak  value  with  AOA  for  a  wing-body-tail  case  than 
for  a  wing-body  or  body-tail  case. 

The  last  physical  phenomena  to  be  modeled  is  the  change  in  axial  force  coefficient  with 
control  deflection  as  AOA  increases.  The  AP95  defines  the  axial  force  term  due  to  control  deflection 
as 


sin  6 


w 


(6A) 


for  the  forward  lifting  surface  and  as 


(Cn. 


W) 


(6B) 
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RGURE  2.  AXIAL  FORCE  AOA  VARIATION  PARAMETERS  FOR 
A  BODY-TAEL  CONHOURATION 


for  a  rearward  lifting  surface  if  two  sets  of  lifting  surfaces  are  present.  If  only  one  set  of  lifting 
surfaces  are  present.  Equation  (6A)  applies  regardless  of  the  fin  location.  Equation  (6B)  contains  the 
tail  interference  term  as  a  result  of  the  forward  set  of  fins,  whereas  Equation  (6A)  does  not  have  this 
term.  Both  and  Cj^  of  Equation  (6)  are  the  normal  force  coefficient  on  the  wing  or  tail  due 
to  both  AOA  an3’ control  3eflection. 

In  analyzing  Equations  (6A)  and  (6B),  they  both  basically  assume  that  the  axial  force  due  to 
control  deflection  is  simply  the  normal  force  of  the  wing  in  conjunction  with  the  body  times  the  sine 
of  the  control  deflection  angle  6.  In  examining  comparisons  of  this  approach  to  data,  it  was  found 
that  when  a  and  6  were  of  the  same  sign.  Equation  (6)  gave  agreement  with  experimental  data  which 
was  quite  acceptable.  However,  when  a  and  6  were  of  opposite  signs,  it  was  found  the  agreement 
was  not  as  good  for  higher  Mach  numbers.  It  is  suspected  that  part  of  the  reason  for  this  is  the 
nonlinear  model  of  and  kB(W)  are  accurate  for  a  and  6  of  the  same  sign,  but  when  a  and  6  are 
of  opposite  signs,  the  total  values  of  kw(B)  and  kB(w)  are  correct,  but  each  may  be  in  error  and  the 
errors  tend  to  compensate.  It  is  speculated  that  the  fundamental  source  of  this  discrepancy  is  wing 
gap  effects  when  the  wings  are  deflected.  This  cancellation  of  errors  could  give  accurate  normal 
force,  but  possibly  inaccurate  axial  force  contributions  due  to  AOA.  The  other  possible  source  of  the 
inaccuracy  in  the  term  is  the  fact  that  kw(B)  does  not  include  nonlinearities  from  control  deflection 
(only  a,  M,  AR,  X  nonlinearities  are  included).  These  nonlinear  wing  gap  effects  could  be  accounted 
for  adequately  when  a  and  6  are  of  the  same  sign  from  the  a  alone  contribution,  but  inadequately 
when  they  are  of  opposite  sign. 
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TABLE  2.  VALUES  OF  THE  AXIAL  FORCE  AOA  PARAMETERS  FOR 
A  BODY-TAIL  CONRCURATION 


M 

f'(M,0) 

f(M,30) 

f(M,60) 

f(M,90) 

0 

0.020 

-0.240 

-0.79 

-0.060 

0.6 

0.020 

-0.200 

-0.66 

-0.060 

0.8 

0.070 

-0.160 

-0.46 

-0.060 

0.9 

0.080 

-0.135 

-0.38 

-0.060 

1.0 

0.120 

-0.110 

-0.26 

-0.060 

1.15 

0.160 

-0.020 

-0.18 

-0.140 

1.2 

0.186 

0 

-0.12 

-0.160 

1.3 

0.200 

0.040 

-0.06 

-0.220 

1.5 

0.258 

0.100 

0.06 

-0.240 

2.0 

0.330 

0.150 

0.20 

-0.130 

2.5 

0.350 

0.180 

0.28 

-0.060 

3.0 

0.346 

0.190 

0.30 

0 

3.5 

0.325 

0.200 

0.32 

0.030 

4.5 

0.230 

0.220 

0.33 

0.065 

2:6.0 

0.180 

0.230 

0.34 

0.070 

Since  no  direct  data  measurements  are  available  and  since  the  control  deflection  matrix  of 
Reference  1  gives  quite  acceptable  values  of  normal  force  and  pitching  moment,  the  approach  taken 
here  is  to  define  a  term  f(M,aw)  to  multiply  Equation  (6)  by  when  a  and  6  are  of  opposite  signs. 

Then 


C.  =  (Cj,,  sin  6^)  f(M,a^) 

^W(B)  ^  ^ 


(7A) 


and 


and  a-j.  of  Equations  (7)  are  defined  by 


(7B) 
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HGURE  3.  AXIAL  FORCE  AOA  VARIATION  PARAMETERS  FOR 
A  WING-BODY-TAIL  CONFIGURATION 


=  a  + 


■'w 


a. 


a  +  6, 


(8) 


Since  the  generic  wind  tunnel  data  bases  available  have  fins  that  are  too  small  to  accurately  determine 
f(M,  a),  use  of  other  data  bases'*’^’  in  conjunction  with  the  AP95  code  will  be  used  instead. 
Figure  4A  and  Table  4A  give  values  of  f(M,aw)  for  Mach  numbers  between  1.5  and  4.60  at  O  =  0, 
and  Figure  4B  and  Table  4B  give  the  complementary  values  at  <&  =  45  deg.  For  Mach  numbers  above 
4.6,  the  values  at  4.6  are  used.  For  Mach  numbers  below  1.5,  f(M,aw)  is  assumed  to  go  to  its  value 
of  1 .0  at  M„  =  0.8  and  remain  at  that  value  below  that  Mach  number.  This  assumption  again  is  based 
on  comparisons  to  data  at  low  speeds.  For  a  and  6  of  the  same  sign,  f(M,a^)  is  always  1.0. 

It  should  also  be  pointed  out  that  the  value  of  f(M,aw)  at  =  90  deg  has  been  assumed  to 
be  0.5  for  both  the  O  =  0  and  $  =  45  deg  planes.  If  the  wing  were  not  in  the  presence  of  the  body, 
then  a  value  of  1.0  would  be  natural  at  least  at  <I>  =0.  However,  due  to  body  interference,  the  full 
effect  of  control  deflection  on  f(M,avv)  at  does  not  appear  to  be  obtained  at  =  90  deg  based 
on  extrapolated  experimental  data  at  lower  values  of  Hence,  the  assumed  value  of  0.5. 
Hopefully,  higher  AOA  experimental  data  will  become  available  that  can  be  used  to  improve  upon 
this  assumption. 
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TABLE  3.  VALUES  OF  THE  AXIAL  FORCE  AOA  PARAMETERS  FOR 
A  WING-BODY-TAIL  CONHGURATION 


M 

f’(M,0)/rad 

f(M,30) 

f(M,60) 

f(M,90) 

0 

0.02 

-0.48 

-0.40 

-0.10 

0.60 

0.02 

-0.40 

-0.54 

-0.15 

0.80 

0.02 

-0.37 

-0.58 

-0.20 

0.90 

0.02 

-0.35 

-0.60 

-0.26 

1.00 

0.03 

-0.33 

-0.61 

-0.38 

1.15 

0.04 

-0.30 

-0.60 

-0.58 

1.30 

0.05 

-0.27 

-0.58 

-0.69 

1.50 

0.09 

-0.23 

-0.55 

-0.72 

2.00 

0.20 

-0.12 

-0.40 

-0.55 

2.50 

0.31 

0 

-0.20 

-0.34 

3.00 

0.42 

0.11 

0 

-0.14 

3.50 

0.53 

0.22 

0.12 

-0.09 

4.50 

0.75 

0.41 

0.23 

-0.09 

^6.00 

0.80 

0.45 

0.26 

-0.09 

Where  Ca(M,90)  =  C^IM.O)  +  f(M,90)  ^  0 


The  final  issue  that  must  be  resolved  is  how  to  break  down  the  changes  in  with  AOA  into 
the  individual  components.  The  AP95'  currently  breaks  these  components  down  into  those  due  to 
axial  pressure  drag,  skin-friction  and  base  drag.  The  present  AP95  methodology  uses  the  work  of 
Reference  13  to  account  for  changes  in  to  AOA  30  deg.  From  AOA  30  to  90  deg,C^  is  assumed 
to  go  linearly  to  zero  from  its  value  at  AOA  30  deg.  The  approach  to  compute  wil!  not  change 
from  Reference  1  since  Reference  1  was  based  on  wind  tunnel  data.  It  will  be  assumed  that  is 
independent  of  AOA.  This  means  that  the  remaining  change  in  with  AOA  is  accounted  for  in  the 
pressure  term.  Thus 


(AC^).  -  f(M.a)  -  (AC^_).  (9) 

Then,  the  individual  components  of  axial  force  of  the  body  are 

“  (^Ap)a=0  ^^^Ap^a  (lOA) 
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TABLE  4A.  f(M,aw)  AT  4>  =  0  deg 


“w 

M  =  1.50 

M  =  2.00 

M  =  2.35 

M  =  2.87 

M  =  3.85 

M  =  4.60 

-20.0 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

-15.0 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

-10.0 

1.0000 

0.9900 

1.0000 

1.0000 

0.9000 

0.7232 

-5.0 

1.0000 

0.9400 

0.8500 

0.7500 

0.6000 

0.2400 

0.0 

0.9200 

0.8200 

0.3000 

0.2500 

0.1400 

0.4000 

5.0 

0.8200 

0.7200 

0.1328 

0.0500 

0.0000 

0.0000 

10.0 

0.7300 

0.6300 

0.1833 

0.0597 

0.0100 

0.0000 

15.0 

0.6700 

0.5700 

0.2800 

0.1100 

0.0200 

0.0000 

20.0 

0.6300 

0.5470 

0.4000 

0.2300 

0.0700 

0.0100 

25.0 

0.5836 

0.5436 

0.4800 

0.3300 

0.1800 

0,0700 

30.0 

0.5803 

0.5403 

0.5000 

0.4300 

0.2700 

0.1600 

35.0 

0.5769 

0.5369 

0.5000 

0.4800 

0.3400 

0.2300 

40.0 

0.5736 

0.5336 

0.5000 

0.5000 

0.4200 

0.3100 

45.0 

0.5702 

0.5302 

0.5000 

0.5000 

0.4700 

0.3600 

50.0 

0.5669 

0.5269 

0.5000 

0.5000 

0.5000 

0.4200 

1 

55.0 

0.5635 

0.5235 

0.5000 

0.5000 

0.5000 

0.4600 

60.0 

0.5602 

0.5202 

0.5000 

0.5000 

0.5000 

0.5000 

65.0 

0.5568 

0.5168 

0.5000 

0.5000 

0.5000 

0.5000 

70.0 

0.5535 

0.5135 

0.5000 

0.5000 

0.5000 

0.5000 

75.0 

0.5401 

0.5101 

0.5000 

0.5000 

0.5000 

0.5000 

80.0 

0.5268 

0.5068 

0.5000 

0.5000 

0.5000 

0.5000 

85.0 

0.5134 

0.5034 

0.5000 

0.5000 

0.5000 

0.5000 

90.0 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

95.0 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

100.0 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 
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TABLE  4B.  f(M,aw)  AT  =  45  deg 


Kw 

M  =  1.50 

M  =  2.00 

M  =  2.35 

M  =  2.87 

M  =  3.85 

M=4.60 

-20.0 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

-15.0 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

-10.0 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

-5.0 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

0.0 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

5.0 

0.9000 

0.9000 

0.9000 

0.9000 

0.9000 

0.9000 

10.0 

0.1000 

0.1400 

0.1800 

0.2200 

0.2600 

0.3000 

15.0 

0.1000 

0.0800 

0.0600 

0.0400 

0.0200 

0.0000 

20.0 

0.1700 

0.1520 

0.1340 

0.1160 

0.0980 

0.0800 

25.0 

0.3100 

0.2740 

0.2380 

0.2020 

0.1660 

0.1300 

30.0 

0.3700 

0.3320 

0.2940 

0.2560 

0.2180 

0.1800 

35.0 

0.4400 

0.3980 

0.3560 

0.3140 

0.2720 

0.2300 

40.0 

0.4800 

0.4400 

0.4000 

0.3600 

0.3200 

0.2800 

45.0 

0.5000 

0.4660 

0.4320 

0.3980 

0.3640 

0.3300 

50.0 

0.5000 

0.4760 

0.4520 

0.4280 

0.4040 

0.3800 

55.0 

0.5000 

0.4860 

0.4720 

0.4580 

0.4440 

0.4300 

60.0 

0.5000 

0.4960 

0.4920 

0.4880 

0.4840 

0.4800 

65.0 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

70.0 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

75.0 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

80.0 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

85.0 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

90.0 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

95.0 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

100.0 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 
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(lOB) 


(CA,)a  -  (CA3)a=0  "  (^C 


(IOC) 


Of  course, 


(CA  =  (Ca,)«  "  (CA^a  -  (CaA  =  C^o  " 


(11) 


It  will  also  be  assumed  that  there  is  no  change  in  the  zero  lift  axial  force  coefficient  of  a  lifting  surface 
with  AOA.  This  assumption  will  also  apply  to  the  components  of  axial  force  on  the  lifting  surface. 
Thus  all  the  axial  force  change  with  AOA  for  bodies  alone  or  those  with  lifting  surfaces  at  6  =  0  will 
be  accounted  for  by  the  body  change  in  with  AOA.  These  changes  are  defined  by  Equations  (9) 
through  (1 1).  The  lifting  surface  change  in  due  to  6  is  accounted  for  by  Equation  (7)  with 
f(M,aw)  being  1.0  when  a  and  6  are  the  same  sign  and  varying  according  to  Table  4  when  of  opposite 
signs. 


3.0  RESULTS  AND  DISCUSSION 


The  first  case  to  compare  the  new  axial  force  methodology  to  experimental  data  is  the  large 
NASA/Tri-Service''*  wind  tunnel  data  base.  Body  alone  wind  tunnel  results  are  compared  to  the  new 
theory  at  Mach  numbers  0.6  to  4.5.  The  wind  tunnel  model  tested  was  12.33  calibers  long  with  a 
3.0  caliber  tangent  ogive  nose  (see  Figure  5A).  It  had  a  boundary  layer  trip  and  was  tested  at 
Reynolds  numbers  of  about  2  x  lOVft  for  most  conditions.  Thus  the  “wind  tunnel  model  with 
boundary  layer  trip”  option  was  used  in  the  aeroprediction  code  to  compute  the  skin-friction 
component  of  axial  force.  This  option  means  that  flow  is  assumed  to  be  turbulent  over  the  entire 
body.  The  wind  tunnel  data  were  given  in  terms  of  forebody  and  chamber  axial  force.  These  two 
values  were  added  together  to  compare  against  the  total  aeroprediction  axial  force  predictions.  It  is 
suspected  that  at  subsonic  Mach  numbers,  the  chamber  axial  force  is  somewhat  different  than  the  base 
drag  with  a  solid  base.  Hopefully,  the  AOA  corrections  will  be  similar  for  the  two,  however. 

Figure  5B  presents  the  results  of  the  comparison  between  the  aeroprediction  axial  force  and 
those  from  wind  tunnel  results  for  the  body  alone  case  of  Figure  5  A.  The  wind  tunnel  data  were  only 
available  to  30  deg  AOA  at  lower  Mach  numbers,  but  up  to  45  deg  AOA  at  high  Mach  numbers.  As 
seen  in  the  figure,  the  new  method  is  superior  to  that  in  the  AP95,  particularly  above  AOA  30  deg. 
It  is  also  seen  that  at  the  higher  Mach  numbers  as  well  as  are  both  predicted  reasonably 
well,  whereas  at  the  lower  Mach  numbers,  C .  has  &e  right  trend  but  total  differs  from  the  data 
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HGURE  5A.  BODY  ALONE  AND  ONE  BODY-TAIL  CONHGURATION  OF 
NASA  TRI-SERVICE  DATA  BASE*'* 
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due  to  .  It  is  believed  the  theoretical  predictions  are  closer  to  actual  fight  conditions  than  the 
model  due°to  the  fact  that  the  model  had  a  large  hollow  chamber  and  the  prediction  assumed  a  solid 
base.  The  base  drag  for  a  configuration  like  this  with  a  solid  base  at  subsonic  Mach  numbers  is  about 
0.12.  This  is  almost  as  large  as  the  total  experimental  axial  force  given  in  Figure  5B  for  M„  =  0.6  and 
0.9.  Hence,  in  viewing  the  comparisons  with  data,  one  should  keep  this  fact  in  mind. 

The  body-tail  configuration  (Figure  5A)  comparisons  of  experiment^"*  and  the  new  method  are 
given  in  Figure  5C  for  the  basic  same  wind  tunnel  parameters  as  for  Figure  5B.  That  is,  a 
configuration  with  a  boundary  layer  trip  tested  at  Rn/A;  of  2  x  10®.  Mach  numbers  shown  in 
Figure  5C  vary  from  0.6  to  4.5.  Figure  5B  utilized  Table  1  and  Figure  1  in  the  new  methodology, 
whereas  Figure  5C  uses  Figure  2  and  Table  2.  As  seen  in  the  various  comparisons  of  as  a  function 
of  AOA  for  various  Mach  numbers,  the  improved  AOA  theoiy  gives  much  better  agreement  to  data 
than  the  AP95  results.  It  is  also  seen  that  the  average  accuracy  goal  of  ±10  percent  on  is  easily 
obtained  with  the  new  theory,  whereas  the  AP95  did  not  achieve  this  goal  at  higher  AOA. 

The  third  case  considered  is  one  of  the  body  alone  configurations  tested  by  Baker.*®  The 
Baker  data  base  was  the  primary  data  base  upon  which  the  theoretical  axial  force  methods  of 
References  3  through  5  were  based.  This  data  base  consisted  of  several  body  alone  and  body-tail 
configurations  tested  at  AOAs  to  180  deg  and  Mach  numbers  0.6  to  3.0.  Figure  6A  shows  the  body 
alone  and  body-tail  considered  as  an  example  here.  The  body  has  a  2.5  caliber  tangent  ogive  nose 
and  7.5  caliber  cylindrical  afterbody.  The  tail  planform  has  an  aspect  ratio  of  one.  The  configuration 
was  tested  with  a  boundary  layer  trip  at  a  Reynolds  number  of  4  x  10‘®/ft. 

Figure  6B  compares  the  new  theoretical  approach  with  the  body  alone  axial  force  data  of 
Reference  15.  Also  shown  on  Figure  6B  for  the  lower  Mach  number  cases  is  the  fourth  order  axial 
force  method  of  Reference  4  and  the  third  order  method  of  Reference  5  for  the  higher  Mach  number 
cases.  Note  that  the  new  method  achieves  one  of  its  objectives  of  being  as  accurate  as  the 
References  4  and  5  methods,  while  using  a  single  versus  dual  methods.  One  of  the  reasons 
References  4  and  5  methods  compare  as  well  as  they  do  to  data  is  the  prediction  that  is  based 
on  the  Baker  data  base.  The  new  method  uses  AP95  theoretical  methodology  to  predict  and  still 
gives  as  good  or  better  comparison  to  the  Baker*®  data  than  either  Reference  3  or  4. 

Figure  6C  compares  the  predictions  of  the  new  theory  to  the  Baker  data  and  the  methods  of 
Ingram®  and  Jorgensen^  at  several  supersonic  Mach  numbers  for  the  body-tail  configuration  of 
Figure  6A.  For  the  Jorgensen  method,  the  AP95  value  of  has  been  assumed.  Note  that  both 
the  third  order®  and  the  present  fourth  order  in  AOA  methods  predict  quite  well  for  supersonic 
Mach  numbers.  However,  the  third  order  AOA  method  does  not  work  well  at  subsonic  Mach 
numbers.  The  Figure  6A  configuration  was  tested  with  a  boundary  layer  trip  and  a  solid  base  so  none 
of  the  chamber  effects  on  base  pressure  are  present  as  was  the  case  with  the  NASA  data  base  of 
Figure  5. 

In  examining  the  comparisons  of  the  new  method  of  this  report  to  those  of  References  3 
through  5  in  Figures  5  and  6,  it  is  seen  that  one  of  the  objectives  of  this  work  has  been  achieved. 
That  is,  a  single  method  (versus  multiple  methods)  has  been  defined  that  works  as  well  as  or  better 
than  those  available®'®  for  body  alone  and  body-tail  configurations  over  the  entire  Mach  number  range 
and  for  AOAs  to  90  deg.  The  method  has  been  validated  for  a  more  limited  Mach  number  range  due 
to  availability  of  data  however. 
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References  3  through  5  derived  nonlinear  AOA  predictors  for  for  body  alone  and  body-tail 
configurations.  No  consideration  was  given  to  configurations  with  more  than  one  set  of  lifting 
surfaces  or  to  control  deflections  so  the  user  of  this  technology  is  left  with  the  use  of  body-tail 
methodology  for  configurations  with  more  than  one  set  of  lifting  surfaces.  The  next  several  examples 
will  apply  the  present  new  method  of  predicting  for  configurations  with  more  than  one  set  of 
lifting  surfaces.  The  theory  for  these  examples  is  base“d  on  AP95  methods  at  AOA  =  0,  Equations  (2) 
and  (3),  along  with  Figure  3  and  Table  3. 

The  first  of  the  configurations  with  two  sets  of  lifting  surfaces  is  a  canard-body-tail 
configuration  with  a  hemispherical  nose  as  shown  in  Figure  7A.  This  configuration  was  tested*’  with 
a  boundary  layer  trip  at  a  Reynolds  number  of  2  x  lOVft.  The  Mach  numbers  considered  were  0.2 
to  4.63  and  AO  As  only  up  to  20  deg.  This  configuration  should  therefore  provide  a  good  test  of  the 
methodology  of  the  AP95  as  well  as  a  check  on  the  initial  trends  of  .  It  is  also  outside  the  data 
base  upon  which  Figure  3  and  Table  3  were  derived. 

Figure  7B  presents  results  of  the  new  method  compared  to  the  data*’  and  the  AP95  for  Mach 
numbers  0.2, 0.6,  1.0, 1.2, 1.75,  2.86,  3.95,  and  4.63.  The  theoretical  computations  are  shown  for 
AOA  to  90  deg  even  though  data  are  given  to  only  20  deg  AOA.  In  comparing  the  AP95  and  new 
method  to  the  data,  it  is  seen  that  near  (a  =  0,  both  methods  give  results  well  within  the  desired 
±10  percent  average  accuracy  goal.  This  is  quite  encouraging  in  view  of  the  fairly  complex 
configuration  of  Figure  7 A.  However,  as  AOA  is  increased  above  0  deg,  the  new  method  gives 
trends  which  in  general  are  superior  to  those  of  AP95,  particularly  at  the  lower  and  higher  Mach 
number  cases  of  Figure  7B. 

The  next  case  considered  is  a  canard  controlled  missile  model  tested  at  the  Naval 
Postgraduate  School  in  Monterey  (See  Figure  8A).  The  configuration  tested  was  a  1/3  scale  model 
which  was  22.6  calibers  in  length.  It  had  canards  with  aspect  ratio  of  1.59  and  fairly  large  tail 
surfaces  with  an  aspect  ratio  of  0.9.  Wind  tunnel  results  are  given  for  AOA  to  45  deg  at  M<„  =  0.2; 
for  both  0  =  0  and  45  deg  roll  positions;  and  for  canard  deflections  of  0,  +20  deg  and  -20  deg. 
Reference  18  stated  that  the  sting  balance  was  based  on  normal  loads  and  therefore  the  axial  loads 
could  have  considerable  error  as  a  result.  The  theoretical  results  were  calculated  based  on  a  boundary 
layer  trip  and  no  boundary  layer  trip  for  control  deflection  of  0  deg.  When  the  forward  controls  were 
defleeted,  it  was  assumed  that  the  flow  would  be  turbulent  over  the  body  surface,  so  the  boundary 
layer  trip  option  in  the  aeroprediction  code  was  used  for  these  calculations.  Freestream  Reynolds 
number  of  1.42  x  lOVft  was  assumed  for  the  computations  based  on  sea  level  conditions. 

Figure  8B  and  8C  compares  the  AP95  and  new  method  to  the  experimental  results  from 
Reference  18.  Reference  18  also  gave  Missile  Datcom*’  results  at  8  =  0,  and  these  results  are  also 
given  on  the  figure.  Figure  8B  gives  the  <I>  =  0  deg  results  and  Figure  8C  gives  the  ^>  =  45  deg 
results.  In  examining  the  ^  =  0  deg  and  6  =  0  deg  results  first,  it  is  seen  that  the  results  for  no 
boundary  layer  trip  appear  to  match  the  data  better  than  with  a  trip  present  and  that  both  the  new 
method  results  (with  or  without  a  trip)  are  superior  to  the  AP95  and  Missile  Datcom  results.  For 
^  =  0  and  6  =  db20  deg,  the  AP95  is  slightly  better  for  the  6  =  +20  deg  case  and  the  new  method  is 
better  for  6  =  -20  deg. 

For  the  $  =  45  deg  roll  and  6  =  0  deg  case  of  Figure  8C,  it  is  seen  that  the  zero  AOA  axial 
force  measurement  agrees  closer  to  the  case  with  a  boundary  layer  trip  than  no  boundary  layer  trip. 
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It  is  also  interesting  to  note  that  this  experimental  value  of  is  0.64  versus  that  in  Figure  8B  of  0.45, 
although  in  principle  the  values  should  be  approximately  the  same  within  measurement  errors.  This 
illustrates  the  point  the  authors  of  Reference  18  made  with  respect  to  the  sting  balance  being  defined 
for  large  normal  forces  resulting  in  inherently  less  accurate  axial  force  measurements.  Note  that  no 
AP95  results  are  shown  on  Figure  8C  due  to  the  fact  the  AP95  was  only  applicable  at  4>  =  0  deg. 
Note  that  the  $  =  45  deg,  6  =  ±20  deg  cases  of  Figure  8C  show  the  new  theory  compares  very  well 
to  experiment.  It  is  particularly  encouraging  to  see  the  theory  predict  the  stall  effect  that  occurs  at 
6  =  +20  and  a  =  20  to  45  deg.  This  is  shown  on  Figure  8C  with  a  fairly  sharp  drop  in  above  AOA 
20  deg. 


A  wing-body-tail  configuration  where  the  wings  are  used  for  control  was  the  next  case 
considered  in  the  validation  process.  This  configuration  is  shown  in  Figure  9A  with  the  experimental 
results  given  in  Reference  20.  This  configuration  has  a  length  of  about  18  calibers  with  a  tangent 
ogive  nose  of  2.25  calibers  in  length.  It  has  wings  and  tails  of  fairly  high  aspect  ratios  of  2.8  and  2.6 
respectively.  Data  were  taken  at  Mach  numbers  1.5  to  4.63,  for  AO  As  to  45  deg  and  control 
deflections  of  0  and  10  deg  at  M  of  1 .5  and  2.0  and  0  to  20  deg  at  M  of  2.35  to  4.63.  The  data  were 
taken  at  a  Reynolds  number  of  2.5  xlOVft  and  boundary  layer  trips  were  also  used.  The  model  had 
a  hollow  chamber,  and  chamber  axial  force  measurements  were  given  separately  in  Reference  20. 
These  results  were  added  to  the  forebody  axial  force  measurements  to  compare  with  the  AP95  and 
new  AOA  axial  force  prediction  method  presented  in  this  report. 

Figures  9B  through  9E  show  the  comparisons  of  the  new  theory  with  the  previous  AP95 
calculations  and  the  data  of  Reference  20.  Mach  numbers  1 .5, 2.35,  3.95  and  4.63  are  shown  in  the 
comparisons.  Figures  9B  and  9C  are  the  0  deg  control  deflection  results  at  O  =  0  and  45  deg, 
respectively.  The  theory  gives  the  same  results  for  these  roll  orientations  at  a  given  AOA  and  Mach 
number.  As  seen  in  the  comparisons  to  data,  there  is  a  slight  effect  of  roll  on  axial  force,  but  not 
enough  to  account  for  in  a  semiempirical  code  for  no  control  deflection.  Figures  9D  and  9E  give  the 
corresponding  comparisons  of  theory  to  data^°  for  <I>  =  0  and  45  deg  roll  respectively,  but  with 
control  deflection.  At  M...  =  1.5,  the  control  deflection  is  10  deg,  whereas  at  the  other  Mach  numbers 
it  is  20  deg.  Here  the  theory  does  distinguish  between  roll  of  0  and  45  deg,  since  only  two  fins  are 
deflected  at  0  deg  roll  whereas  all  four  are  deflected  at  45  deg  roll.  Also  note  that  no  AP95 
computations  are  shown  on  either  Figure  9C  or  9E  due  to  the  fact  the  AP95  is  only  applicable  at 
0  =  0  deg.  Note  that  in  all  the  Figures  9B  through  9E,  the  new  method  in  general  is  equal  to  or 
superior  to  the  AP95  and  also  gives  very  good  agreement  in  comparison  to  data.  The  worst 
comparison  is  for  M  =  3.95,  6  =  20  deg,  0  =  0  deg  and  AOA  40  deg.  Even  for  this  worst  case 
comparison,  the  new  theory  differs  from  data  by  only  12  percent.  The  average  accuracy  errors  of  the 
new  method  to  data  in  Figures  9B  through  9E  are  well  within  the  ±10  percent  accuracy  level  claimed 
in  the  overall  accuracy  levels  for  normal  and  axial  force  coefficients. 

The  final  configuration  considered  for  validation  is  a  wing-body-tail  case  where  the  tail 
surfaces  are  used  for  control.  The  configuration  is  identical  to  configuration  9A  except  the  wing 
trailing  edges  are  tmncated  with  a  height  of  0.05  inches  versus  being  sharp,  as  shown  in  Figure  9A. 
This  means  the  drag  of  this  configuration  is  slightly  greater  than  the  wing-controlled  case  due  to  the 
wing  trailing  edge  base  pressure  drag.  The  wind  tunnel  data^'  were  taken  at  a  Reynolds  number  of 
2.0  X  lOVft  versus  the  wing  controlled  case  of  2.5  x  lOVft.  However,  a  boundary  layer  trip  was  used 
for  both  the  tail  and  wing  controlled  models. 
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HGURE  9A.  AIR-TO-AIR  MISSILE  CONHGURATION  USED  IN  VAUDATION^ 
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nOURE  9B.  COMPARISON  OF  THEORY  AND  EXPERIMENT  FOR  CONHGURATION  OF 

HGURE  9 A  ($  =  0  deg,  6  =  0  deg) 
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HGURE  9D.  COMPARISON  OF  THEORY  AND  EXPERIMENT  FOR  CONFIGURATION 

OF  HGURE  9A  (<I>  =  0  deg) 
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HGURE  9E.  COMPARISON  OF  THEORY  AND  EXPERIMENT  FOR  CONHGURATION 

OF  HGURE  9A  (4>  =  45  deg) 
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The  comparisons  of  axial  force  coefficients  to  data^'  of  the  new  method  to  the  AP95  and 
experiment  for  no  control  deflection  is  quite  similar  to  the  results  of  Figures  9B  and  9C  and  therefore 
will  not  be  repeated  here.  The  results  for  6  =  -20  deg  at  <I>  =  0  and  €>  =  45  deg  are  shown  in 
Figures  lOA  and  lOB  respectively.  Results  are  given  in  each  figure  for  Mach  numbers  of  1.  5, 2.87, 
3.85,  and  4.6.  The  theory  is  shown  up  to  AOA  90  deg  although  data  are  available  only  up  to  AOA 
of  35  deg.  For  Mach  number  of  1.5,  both  the  AP95  and  the  new  method  appear  to  give  acceptable 
results  compared  to  data  at  both  0  =  0  and  O  =  45  deg.  However,  for  the  three  higher  Mach 
numbers,  the  new  method  is  superior  to  the  AP95  at  both  roll  positions.  The  reason  for  this  improved 
performance  is  due  to  the  empirical  model  for  AOA  and  control  deflection  of  opposite  signs  as  given 
in  Figure  4. 


4.0  SUMMARY  AND  RECOMMENDATIONS 


An  improved  semiempirical  method  for  axial  force  calculation  on  missile  configurations  has 
been  developed.  The  method  uses  the  theoretical  methods  currently  used  in  the  1995  version  of  the 
NSWCDD  Aeroprediction  Code  (AP95)  for  zero  AOA  axial  force  calculations  and  wind  tunnel  data 
bases  to  compute  changes  in  axial  force  at  AOA  for  body  alone  and  wing-body  configurations. 
Comparisons  to  data  and  existing  theoretical  approaches  for  computing  the  AOA  contribution  to  axial 
force  indicate  the  new  method  to  be  as  good  as  or  superior  to  existing  techniques  at  all  Mach  numbers 
and  AO  As.  The  improved  method  was  also  extended  to  configurations  which  have  two  sets  of  lifting 
surfaces  and  also  to  consider  axial  force  due  to  control  deflection.  Comparisons  of  the  extended 
method  to  the  AP95  for  wing-body-tail  configurations  showed  improvements  to  the  AP95  at  higher 
AOA,  higher  Mach  number  and  when  the  AOA  and  control  deflection  were  of  opposite  sign. 

In  developing  the  improved  semiempirical  method  several  assumptions  were  necessary.  These 
assumptions  include: 

a.  The  value  of  C^  at  AOA  90  deg  for  configurations  with  2  sets  of  lifting  surfaces  and  for 
configurations  with  control  deflections. 

b.  The  value  of  C^  at  high  AOA  outside  the  Mach  number  range  for  which  data  were 
available  (M  >  4.6). 

c.  How  to  break  the  contribution  in  C^  due  to  AOA  into  individual  axial  force  components 
due  to  skin-friction,  wave  or  pressure  and  base  pressure. 

As  a  result  of  these  assumptions,  it  is  recommended  that  Navier  Stokes  (N.S.)  computations 
be  conducted  on  a  limited  set  of  cases  to  try  to  sort  out  the  axial  force  components.  These 
computations  should  be  conducted  in  conjunction  with  appropriate  wind  tunnel  tests  at  high  AOA 
for  verification  of  the  N.S.  computations.  It  is  then  recommended  that  the  current  methodology  be 
refined  based  on  results  from  these  wind  tunnel  data  and  N.S.  computations. 
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FIGURE  lOA.  COMPARISON  OF  AXIAL  FORCE  COEFFICIENTS  OF  A  TAIL  CONTROLLED 
WING-BODY-TAIL  CONFIGURATION  (0  =  0  deg) 
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6.0  SYMBOLS  AND  DEFINITIONS 

Angle  of  Attack 

1972  and  1995  versions,  respectively,  of  the  NSWCDD 
aeroprediction  code 

Axial  force  coefficient 

That  part  of  the  axial  force  coefficient  due  to  the  base  of  the 
configuration 

That  part  of  the  axial  force  coefficient  due  to  skin  friction 

Axial  force  coefficient  at  AOA  of  zero  degrees 

That  part  of  the  axial  force  coefficient  due  to  pressure  on  the  body 
excluding  the  base  region 

Axial  force  coefficient  component  due  to  AOA 

Axial  force  coefficient  component  due  to  control  deflection,  6 

Drag  coefficient 

Lift  coefficient 

Function  and  its  derivative  respectively,  which  are  defined  based  on 
various  data  bases  as  a  function  of  Mach  number  and  AOA 

Lift  to  drag  ratio 

Mach  number  and  freestream  Mach  numbers 
Modified  Newtonian  Theory 
Angle  of  Attack  (degrees) 

Control  deflection  (degrees)  with  leading  edge  up  as  positive 

Roll  position  with  0  =  0  deg  having  fins  in  plus  (+)  roll  orientation 
and  O  =  45  deg  having  fins  in  cross  (x)  roll  position 
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